Effects of relaxation processes during deposition of anisotropic grains on a flat 

substrate 
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The ballistic deposition on a one dimensional substrate of grains with one degree of freedom, called 
spin, is studied with respect to relaxation processes during deposition. The "spin" represents the 
grain anisotropy, e.g. its longest axis with respect to the vertical. The grains interact through some 
contact energy (J) and are allowed to flip with a probability q during deposition and relaxation. 
Different relaxation processes are investigated. The pile structure is investigated, i.e. the density 
and "magnetisation", as a function of q and J. A percolation transition is found across which the 
cluster size changes from exponential-like to a power law-like dependence. The differences between 
"ferromagnetic" and "anti- ferromagnetic" -like contact energies are emphasized as a function of q. 



I. INTRODUCTION 

The sandpile problem Q] is one of the most recently 
often tackled problems in the self-organisation of com- 
plex systems. There are several questions of interest, not 
only concerning the density, compaction, stability, intrin- 
sic dynamics, etc. but also concerning the detailed non 
equilibrium physics inherent to such systems 0. Exper- 
imental observations pertain to the force structure [3j in 
the pile and to the clustering or percolation properties of 
rain piles Q . It is known that the system is hyperstatic 
In fact the systems are not at all in equilibrium in a 
thermodynamic sense, resulting in hysteresis and ageing 
processes. A description of the construction of the pile is 
thus a relevant input before further studies. 

One important physical constraint to be considered 
in describing granular piles is that the materials are 
not made of symmetrical (spherical or cubic) entities^. 
Whence grains can be imagined to be identical entities 
but having one degree of freedom, call it a spin, for usual 
statistical mechanics considerations. Such a "spin" al- 
lows us for referring to a direction or a rotation process. 
The spin indicates, e.g. the orientation of the longest 
grain axis with respect to the vertical. 

We have studied such a model 0, Q , within a (mag- 
netic) ballistic deposition process; it is similar to the 
Tetris model [8|, but more simplified. Basically the (2- 
dimensional) grains are anisotropic and are deposited 
vertically as in a gravity controlled rain. For simula- 
tion studies we have imagined a two dimensional vertical 
silo like a triangular lattice with vertical edges, but with 
periodic boundary conditions. The grains can be submit- 
ted to an external field during deposition such that their 
orientation can change (or not) at each step, with a cer- 
tain probability. The grain energy is calculated at each 
falling down step according to a sort of Ising energy, the 
exchange integral being a map of the contact energy be- 
tween grains. The latter can have various origins: surface 
roughness, local charges, mechanical, chemical, dipolar 
or magnetic effects, ... . It is known that (long range) 



magnetic interactions lead to clustering and aggregation 
0. For simulation ease, we can limit the interactions to 
a short range. 

We have previously examined clustering effects as in 0, 
through the cluster fractal dimension, the pile overall 
density, the distribution of grain orientations, ... and 
have found through percolation arguments 0, [n| the 
existence of different cluster growth " mechanisms" . 

However, except for a Metropolis like dynamical con- 
dition during the Monte Carlo simulation, we have 
not fully let the system relax during the deposition. This 
constraint is hereby removed. We have performed simula- 
tions under two pile relaxation conditions to be described 
in Sect. [H] One is a total gravitationally controlled re- 
laxation, the other is a blocked gravitationally controlled 
relaxation, relaxing the ballistic constraint, and allowing 
for the grain to move to a neighboring column. It should 
be expected that such a relaxation has a coarsening effect 

The results are found in Sect. IIIII while Sect. IIVI con- 
tains some brief conclusion 



II. 



EXPERIMENTAL PROCEDURE 



A few changes to the experimental procedure presented 
in ||j are hereby introduced. We propose an algorithm 
for creating a pile with a fixed probability q £ [0; 1] for 
spin flip. The algorithm for arbitrary q (so called q-MBD 
model, in contrast to the i-MBD model @) goes as fol- 
lows: 



1. first, we choose a 1-dimensional crenellated hor- 
izontal substrate of spins with a predetermined 
(hereby random) configuration; periodic boundary 
conditions are used; on each site a vertical column 
of empty cells is placed. The lattice symmetry is 
thus hexagonal. 

2. a falling (up or down) spin is dropped along one of 
the lattice columns from a height r max + 5a, where 
fmax is the largest distance between an occupied 
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cluster site and the substrate, i.e. the height of the 
highest column growing from the substrate on the 
lattice 

3. more importantly at each step down the spin can 
flip i.e., change its "sign", with probability q; more 
specifically the probability to choose the "up" di- 
rection is q; 

4. the spin goes down flipping until it reaches a site 
perimeter of the cluster at which time the local gain 
in the Ising energy (in absence of an external field) 

f3E = -(3 J (!) 

<i,j> 

thus for possible cluster growth is calculated. The 
exchange integral J describing spin-spin interaction 
can be thought in granular matter to be like some 
contact energy due to the surface roughness or the 
clastic, electric, magnetic, chemical ... character 
of grains; it can be also a measure of some mis- 
orientation entropy. If the energy gain is negative 
the spin sticks to the cluster immediately (sticking 
probability =1.0) and one goes back to step J3J. In 
the opposite case the spin sticks to the cluster with 
a rate exp(— A/3E) where A/3E is the local gain 
in the Ising energy. If the spin does not stick to 
the cluster it continues going down under the same 
stick and fall rules. At some step the spin sticks to 
the cluster. If the site just below the falling spin 
is occupied the spin is allowed to relax by moving 
down in one of the nearest columns. 

5. At this point we introduce and study three kinds 
of relaxation processes: 

I. No relaxation - the no-relaxation process is 
the pure g-MBD, 0, i.e. the spin sticks to 
the cluster immediately. 

II. Total gravitationally controlled relaxation - 
in the total gravitational relaxation process, 
when the spin reaches the cluster, i.e. the site 
just below the spin is occupied, the spin can, 
with equal probability: stay on the site, or 
roll down all the way in the left or the right 
neighbouring column if one of the first neigh- 
boring sites is unoccupied. The spin continues 
rolling down until it reaches a site from which 
the rolling is not possible; one assumes that 
the rolling process is due to the gravity with- 
out any resistance. 

III. Blocked gravitationally controlled relaxation - 
after reaching the cluster the spin can, with 
equal probability: stay on the site, roll down 
on the left or on the right from the site as for 
the total gravitational relaxation (II) , but can 
also stick to the pile on its way to the lowest 



level: During every move towards the lowest 
level the falling spin can be stopped and stick 
to the cluster depending on the local gain in 
the Ising energy (|TJ. 

If the spin cannot move (due to the relaxation), it 
immediately stops and sticks to the cluster. 

6. After the spin sticks to the cluster one goes back to 
step ©. 

After dropping a (large) number of spins the physi- 
cal quantities of interest like magnetization, density, his- 
togram of up/down spin cluster size, histogram of hole 
cluster size, etc. are computed. 

III. NUMERICAL RESULTS 

All results reported below are for a triangular lattice 
with vertical edges of horizontal size L = 100, i.e. the 
width of the seed substrate, and when the pile made of 
clusters has reached a 100 lattice unit height. The sub- 
strate consists of spins with either direction chosen with 
probability q. Every reported data point hereby corre- 
sponds to an average over 1000 simulations. 

A. Density 

We define the density of a cluster as 

number of spins in the cluster 
^ number of sites on the lattice ' 

in which obviously the number of lattice sites in the de- 
nominator = 10000. 

Figs DEI illustrate the behaviour of the density with re- 
spect to the q and (3J parameters. The figures convince 
us that the density reached by the pile in the presence of 
relaxation is always higher than in absence of relaxation, 
at fixed q and f3J values, as expected. The reached den- 
sity values are not trivial though (compare with 0). It 
can be noticed, as expected also that p n > p nl > p 1 
whatever the interaction sign. 

B. Magnetization 

The dependence of the magnetization defined as 



n+ + n_ 

is shown in Fig|3Ka) as a function of (3 J and q, where 
n + and rt_ are the number of up and down spins re- 
spectively, i.e. n + = 10000,0+. This quantity can be 
considered as a measure of the difference in grain orien- 
tations in the packing, when the spin (+) is understood 
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FIG. 1: Projections of the total spin (or grain) density dependence on q and f3J where flj > 0: (a) case I, (b) case II, (c) case 
III. 
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FIG. 2: Projections of the total spin (or grain) density dependence on q and j3J where /3J > 0: (a) case I, (b) case II, (c) case 
III. 



e.g. as indicating that the grain longest axis orientation 
is along the vertical. 

In [|| a non linear-like dependence on q was found for 
the magnetization though with small departure from lin- 
earity due to the changes in the interaction strength. No- 
tice that the relaxation processes make the departures 
from linearity still smaller; the smallest departures exist 
in the case II, total gravitational relaxation (FigO)) . As 
it was mentioned in [f| this has to be expected if one 
recalls that q, in a first approximation, is the fraction of 
deposited up-spins. 



C. Percolation 

In this section we present results concerning the perco- 
lating cluster (from the substrate to the top of the pile) 
of a typical entity, i.e. the up-spins. In order to do so ev- 
ery pile built by the algorithm has been checked in order 
to determine whether there is a percolating cluster of up 
spins. At fixed q and j3J the fraction of piles consisting 
of this given type of percolating cluster was computed. 
This value is thought to be representative of the critical 
percolation p c . 

Fig0| shows the behaviour of p c with respect to q for 
several values of the j3J parameter for the three cases. 
The sigmoid curve is likely due to the finite size of the 
system. Observe that the interval of q for which < 
p c < 1 is wider when relaxation occurs. One can also 
observe that the cases AF or F and f3J = are much more 
distinguishable when there is relaxation (total or blocked 



gravity). However some convergence toward some finite 
p c value seems to occur for increasing \@J\, as seen on 
FigEl for case II, - the asymptotic value depending in 
each case on (3 J sign and being at a lower q value when 
pj > 0. 

Notice that for (3J > in the relaxation cases per- 
colation occurs for q > 0.59 = qo (Fig^Jj) while in the 
non-relaxation case it exists for q > 0.8 = qo (FigEt). 
Therefore, relaxation has an obvious essential influence 
on the critical percolation value. Recall that in a first 
approximation the q parameter is the percentage of up 
spins in the pile, therefore a lower value of qo means that 
one does not have to use a lot of up spins to have a per- 
colation path (this is useful for example if the up spin 
material is much more expensive or dangerous). 

D. The size (mass) distribution of clusters 

Many problems in granular matter concern the cre- 
ation of clusters according to se greg ation p ^ . creation of 
avalanches 1 15j or decompactionflq. In this section wc 
present the results of the analysis of the size (mass) of 
spin and hole clusters in this q-MBD model with or with- 
out relaxation processes. 

In it has been found that the number (N s ) of clus- 
ters dependence on the size (s) of the clusters demon- 
strates a crossover effect, i.e. from an exponential-like to 
a power law-like dependence. The power law-like regime 
exists for high q values and the exponential law regime 
occurs when the q parameter is low. It was conjectured 
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FIG. 3: The magnetization for (a) the relaxation case I, (b) the relaxation case II, (c) the relaxation case III. Three values of 
f3J are presented. 




FIG. 4: Dependence of the percolation probability p c on q: (a) the relaxation case I, (b) the relaxation case II, (c) the relaxation 
case III. Three values of ft J are presented. 
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FIG. 5: Dependence of the percolation probability p c on q for the relaxation case II (total gravitationally controlled relaxation) 
for: (a) /3J > and (b) (3J < 0. 
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FIG. 6: The exponent k p form power law regime: (a) the relaxation case I, (b) the relaxation case II, (c) the relaxation case 
III. 
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that the border value is at the critical percolation (q c ). 

A crossover effect likewise occurs in presence of relax- 
ation. The whole data analysis is not reproduced here; 
it was done as for the case I in 0. The differences are 
particularly noticeable in the power law region i.e. for 
high values of q. FigH3 shows the behaviour of the power 
law exponent k p , where N$ oc e~ k v s and N s is a number 
of clusters with mass s for the three models. Notice that 
the k p parameter is quasi the same (taking into consider- 
ation error bars) for finite (3 J and (3 J = 0, when q > 0.8, 
i.e. whenp c w 1. In that CELS6 kp is almost q independent. 

Let us introduce (like in 0) two parameters to charac- 
terize the exponential law regime: an exponent &e where 
N(s) ~ exp(— fcss) and Ae where &e = exp(— Ae<z)- Ob- 
serve (FigEJ in contradistinction with the no relaxation 
models that Ae ^ const for (3 J < 0, but Ae — const for 
(3J < —2. In other words Ae = const for \j3J\ < 2 in 
the presence of relaxation, i.e. for small values of q the 
cluster growth does not much depend on the interaction 
strength for higher interaction vaules. 

E. The analysis of spin configurations 

The probability of spin sticking for a given spin neigh- 
bourhood configuration has been counted. The results 
reveal that the most probable configurations depend on 
q and (3 J. We have investigated many configurations in 
the (q, (3 J) plane. We confirm that this approach allows 
to estimate as well the percolation transition, but such 
results cannot be fully displayed here. 

The probabilities seem to have a very small dependence 
on the f3J parameter. However the q parameter induces 
a strong influence on the configuration probability occur- 
rence. No simple empirical fit has been found for such 
dependences. For the sake of conciseness and argument 
we only display two tables (Table III and HT|) for (3, J = — 5 
and for j3J = 5 respectively for case II and various q 
values. 

We notice that the same configurations are most prob- 
able, in both relaxed and non-relaxed cases. The prob- 
ability of occurrence for the most frequent configuration 
when (3 J < is slightly smaller for the non-relaxed case, 
but this implies that the other configurations have a 
higher probabilities and occur more often when there is 
no relaxation. This entails a larger diversity of possi- 
ble configurations in absence of relaxation proceses. For 
(3J > the probability of occurrence for the most fre- 
quent configuration in the non-relaxed case is not neces- 
sarily smaller compared to the II and III case; the dif- 



ferences are slender and may be sensitive to numerical 
details. 



IV. CONCLUSIONS 

Granular pile spreading processes are driven by coop- 
erative non-linear evolution rules. This leads to develop 
patterns which often reach a high level of complexity. 
We have combined topology and mass (or size) in or- 
der to describe granular piling in a simple way on a 1-D 
substrate. The grains are taken to be anisotropic, and 
the pile forming through a ballistic deposition process 
at first. However, in order to take into considerations 
and observations by wonderful colleagues from the phar- 
macy department [13, we have allowed the granular pile 
to somewhat relax, in removing the constraint of ballis- 
tic deposition when the grain reaches the pile. We have 
allowed for its downward motion along neighbourhood 
columns. 

The pile density of course increases with relaxation. 
The size distribution of clusters follows intricate laws. 
Above a percolation transition, a power law distribu- 
tion is found, but it is like an exponential law below. 
Such a transition occurs at non trivial values in the (q, 
(3J) plane. Slight cluster configuration differences exist 
whether the "spin-spin interaction energy" is positive or 
negative. 

The above changes in density could be usefully cross- 
checked with new experimental studies. This would give 
some information on the constact energy (3 J and the spin- 
spin probability, the latter through the measurment on 
the falling spin. Time dependent effects should now be 
usefully investigated both in simulations and in exper- 
imentally realistic cases. For example Majumdar and 
Dean 18] found that the magnetization relaxes extremely 
slowly, in an inverse logarithmic fashion for 1-D chains. 
The higher dimesnionality cases should be of interest. 
The role of avalanches, as in real sand/grain piling should 
also be of interest next, even taking into account the top- 
pling of clusters of granules [HSf. 
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FIG. 7: The exponent Ae form exponential law regime: (a) the relaxation case I, (b) the relaxation case II, (c) the relaxation 
case III. 



TABLE I: The distribution of typical configurations in the case (II) for /3J = 5; con], type corresponds to a configuration type 
(see Fig|HJ; prob. corresponds to a probability magnitude for the configuration type. The symmetrical (with respect to vertical 
axis) configurations are not recognized, i.e. 0-000 and 000-0 are the same configuration types. 
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TABLE II: The distribution of typical configurations in the case (II) for {3J = —5; conf. type corresponds to a configuration 
type (see FiglHJ; prob. corresponds to a probability magnitude for the configuration type. The symmetrical (with respect to 
vertical axis) configurations are not recognized, i.e. 0-000 and 000-0 are the same configuration types. 
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FIG. 8: The label notations for Tables Q and El the arrows denote a spin sign (+ is up), if on the site there is any spin; an 
empty site is "0", the black circle denotes the sticking spin. On the right for both pictures the corresponding configuration 
type is written. 



